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Abstract 

We derive the bulk viscous time scale of neutron stars with quark 
matter core, i.e. hybrid stars. The r-mode instability windows of the 
stars show the theoretical result accords with the rapid rotation pulsar 
data. The fit gives a strong indication for the existence of quark matter 
in the interior of neutron stars. Hybrid stars instead of neutron or strange 
stars may result in submillisecond pulsars if they exist. 

PACS numbers: 97.60. Jd, 12.38.Mh, 97.60. Gb 

Searching for quark-gluon plasma(QGP) and quark matter stars are of great 
interest in high energy physics as well as astrophysics. Compact objects are 
thought to provide a unique astrophysical environment for the transition from 
hadron matter into quark matter. A neutron star with quark matter core is 
usually called hybrid star. However the prediction of distinguishable signature 
from neutron stars is urgent and difficult task because of similar in mass and 
size for the two kinds of stars and their same surface properties as well. An 
alternative way is to clarify limiting rotation constrained from the dynamics 
of dense matter in the interior of the stars. Recently the discovery of r-mode 
instability in relativistic stars jjj gives us a chance to realize this issue. A series 
of papers have investigated the many implications for gravitational radiation 
detection and the evolution of pulsar §, || ||, |, |, |, |, [HJ fy], [[§. Two 
crucial issues regarding the star evolutionary relevance of the r-mode instability 
also have been solved: key results concern the interaction between oscillations 
in core fluid and the crust (l3[ |l5| and bulk viscosity coupled to r-modesfl6f. 

We have known two evident characteristics regarding rotation of pulsars 
from the observed data, the most rapidly spinning pulsar with 1.5ms period 
and the clustering of spinning frequencies of the low mass x-ray accretion bina- 
ries(LMXBs). Many investigators have made efforts to recognize the observable 
properties or what they indicate in neutron star or strange star models through 
the studies on so called r-mode instability window. In the beginnings, some 
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works find hot, young neutron stars would be spun down to 15ms-20ms from 
their birth due to gravitational wave radiation because of r-mode unstable^. 
Soon after, neutron stars with solid crust could be thought a candidate model 
to explain the rapidly rotating accreting pulsar's data through thermal run- 
away recycle of r-mode when the viscous boundary layer damping taken account 



into 1 13 UM. Meanwhile, it was proposed that newly born strange stars would 
also be spun down to about 2.5ms-3ms due to a huge bulk viscosity relative to 
neutron star matter 

However, the resulting viscous heating in the boundary layer for neutron 
stars with solid crust is so intense that it can heat the crust-core interface to 
the melting temperature of the solid crust if the r-mode amplitude is larger than 
some critical value, which is crudely estimated to be ~ 10~ 3 by OwenJijJ and 
a more accurate value is given as 5 x 10~ 3 by Lindblom et al.|Q. Therefore 
Lindblom et al. found it appear unlikely in neutron stars that the r-mode 
instability is responsible for limiting the spin periods of LMXBs]l8||. For strange 

stars, the spin-down to about 3ms corresponds to temperature in 10 6 10 7 K 

instead of the inferred core temperature as a few times 10 8 K, which is a serious 
drawback. It isn't genuinely convincing if the drawback is simply attributed to 



under-representing statistics due to a low number of objects 11 . 

In this paper, we will focus on the r-mode instability window of neutron 
stars with quark matter core and investigate the spin-evolution of the stars in 
presence of r-mode instability . As demonstrated below, the rapidly rotating 
pulsars would be a strong indication of the existence of quark matter in the 
interior of neutron stars(i.e., hybrid stars). A hybrid star would also be inferred 
to a submillisecond pulsar if the quark matter core is large enough. 

In rotating relativistic stars, gravitational radiation drives the r-mode while 
various dissipation mechanisms counteract the fluid motion. In general, shear 
viscosity, or surface rubbing if the star has a solid crust, suppresses the mode 
at low temperatures but bulk viscosity dominates at high temperatures. The 
critical rotation frequency for a given stellar model as a function of temperature 
follows from 

_L + _L + _L + _L = , (1) 

Tgr 7~ sv Tfr v T sr 

where r gr < is characteristic time scale for energy loss due to gravity wave 
emission, t sv and Tb v are the damping times due to shear and bulk viscosities, 
r sr is surface rubbing due to the presence of the viscous boundary layer. 

A polytropic equation of state with a low index n is a good approximation 
for compact stars, as was discussed in many papers^, |l6). The time scale for 
n = 1 polytrope from gravity wave emission can be written 

r gr = -3.26{irGp/n 2 ) 3 (2) 

where G is the gravitational constant, Q is angular rotation frequency, and p 
is the mean density. Since both the shear viscous damping and the rubbing 
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dissipation due to viscous boundary layer dominate at low temperatures and 
some investigations have shown t sv 3> r sr ]l3|, Q , only the surface rubbing need 
to be considered. In refjli], the time scale is estimated as 

r - - 200 x (ok) 2 (t^c) 2 (3) 

We now need to evaluate the bulk viscous damping time scale of hybrid 
stars. We can follow to work out bulk viscosity coupled to r-modes taking 
respectively the viscosity coefficients of strange quark matter(SQM) and neutron 
star matter in the core and the outer covering of the star, but the calculation is 
a complicated and arduous task. In fact, we will only consider the dissipation on 
r-modes in quark matter core because of a huge bulk viscosity of SQM relative to 
neutron star matter. As demonstrated below, we can here avoid the complicated 
calculations to deduce the time scale with a simple method following the formula 
in 0. 

For the given energy of modes E and its variable rate E = dE/dt, the time 
scale is defined as 

- = (4) 
t IE K ' 

The r-mode energy in a rotation relativistic star with radius R has been ex- 
pressed as 



E = ■„ 
2m 



{m + l) 3 (2m + l)\R 4 n 2 J p{r) (-J dr + <3(ft 4 ), (5) 



where a is the dimensionless amplitude, m is natural number and p(r) is the 
density function. The energy dissipation rate due to bulk viscosity can be 
evaluated with the following formula 

E bv = - I (5cr6cr*d 3 x, (6) 



where C is bulk viscosity coefficient, 5a is deformed velocity tensor. We first 
decompose the energy E into two parts, i.e., the quark core and nuclei envelope 
denoted respectively by c and e 

/ T> \ 2m-2 

E ={-£) (Ec + E e ) (7) 

with 

2 rRa / \ 2m+2 



E c = —(m + iy(2m+l)\Rtn'j^ p c {r) I — ) dr, (Si 



2 pR / \ 2m+2 

E e = -^(m+ If {2m + l)!^fi 2 / p e (r) — dr. 



(9) 
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Assuming uniform core and envelope that is approximately correct for neutron 
stars, we immediately estimate 



Pc \Rc 



2 m +3 



Thus the energy can reads 
E = 



Pe ( B_ 

Rr 



2m+3 



Er. 



(10) 



(11) 



An analogous treatment can be made for the energy dissipation, but we take 
the estimation, da ~ (tj-)™ R 2 c (if)™ l+ > into consideration according ]l6| . 
Regarding the fact that the SQM bulk viscosity dominates [^0|, i.e., Csqm 
Cnsm, we easily obtain the dissipation due to bulk viscosity 



Therefore we derive the time scale as 



E Bc 



M -M r f R 



Mr 



2ni 



TBc, 



(12) 



(13) 



where M c represents the mass of SQM core in the interior of stars. Evidently, the 
time scale tb c is in agreement with that of a strange star with radius R Cl which 
can be calculated @. For given R, R c , M, M c (10km, 5km, 1.4M©, 0.7M Q , we 
obtain 



Tbv 



5.4 x 10^s 



M 



IAM Q V 10km / 



(—)' 



(- 



m s 



V lOOMcV 



1000 



(14) 

where m s is strange quark mass. In fact, the reasonable radius of quark matter 
core in interior of hybrid stars is able to range from 3km to 6km The time 
scale strongly depend on the radius of the quark matter core from Eq(13),it is in 
the range of 2.4 x 10 3 s - 4.7 x 10 1 s under the situation of the scaling quantities 
given in the above formula. 

To have those time scale quantities, we can apply equation(l) to evaluate 
the critical frequency of the rotating stars as a function of temperature. Figure 
1 plots the instability window for the given parameters: R — 10km, R c = 5km, 
M = 1AM Q ,M C = 0.7Af Q ,m s = 200MeV. The result(solid line) in figure 1 
shows that the remain pulsars in supernova explosions will reach the instability 
window due to cooling and then slow down by gravitational wave emission brak- 
ing (high braking index 9) along the right edge of the window. This is similar 
to bar strange stars but the stars here only spin down to 704Hz (a period of 
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1.42ms). Especially, the lowest position of the window can be located at the 
inferred star core temperature, truly consistent with the data unlike situation 
of bar strange stars. In contrast to neutron star (without quark matter core), 
our new results show the pulsars follow a track indistinguishable from the crit- 
ical rotating curve govern by the competition between bulk viscous dissipation 
and gravitational wave emission to be spun down according to the mechanism 
discussed by Andersson et al.Q and MadsenQ. Obviously, the solid crust 
in hybrid stars will never melt even if the stars accreting spin-up. To display 
clearly hybrid star advantage, we have also collected the results for neutron 
stars and bar strange stars in figure 1. Finally, it is worthwhile discussing the 
dependence of limiting rotation of pulsars on the size of quark matter core. 
Figure 2 shows the limiting rotation frequency of pulsars as hybrid stars in the 
range of 637-826Hz in temperature region of 1 x 10 8 — 4 x 10 8 K, based on sup- 
position that the radius of the quark matter core takes 3-6km, as seen above, 
from nuclear theoretical prediction The distributed data region are full in 
agreement with our those results. 

In addition, we can imagine that the existence of the higher rotation pulsars 
is possible if there would exist the neutron stars with larger quark matter core. 
The thick line in figure 1 indicates the upper limit of the rotation of the pulsars 
as hybrid stars, which is higher than neutron stars and strange stars and over 
kHz (a period of submillisecond). This implies that submillisecond pulsars, if 
they exist, should be hybrid stars instead of neutron stars or strange stars, which 
needs future study. 

In conclusion, the bulk viscous time scale for hybrid stars has been calculated 
for a range of parameters. The r-mode instability window has been obtained. 
The rapid rotation pulsar data are a significant indication of the existence of 
hybrid stars instead of neutron stars or strange stars. Of course, this is a pos- 
sible evidence of the existence of Quark Matter in the interior of neutron stars. 
Meanwhile, hybrid star model may predict the existence of submillisecond. 
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